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Abstract. In this paper, we characterize bounded ancient solutions to the time-dependent 
Stokes system with zero boundary value in various domains, including the half-space. 

1. Introduction 

In this paper, we show that any solution u(x,t) G L°°(Q x (— oo, 0)) to 

d t u — Au + Vp = 
div u = 



in Vt x (— oo, 0), and (1.1) 
u\ 9Q = (1.2) 



satisfies 

if Q C R n is a bounded domain and n > 2, 

u(t) if Q = R n and n > 2, 

• u(t, x n ) with u n = if O is half space x n > and n > 2, \ • ) 

a(t) + 0( | g |n- 2 ) as |x| — >■ oo if f2 c i? n is an exterior domain with n > 3. 

Throughout the paper we assume that the domains Q have smooth boundary. One has to 
be somewhat careful with the definition of the boundary condition u\gn = since a-priori we 
only assume u to be bounded, with no further regularity assumptions. The usual definition 
is the following: 

Definition 1.1. We call u(x, t) G L°°(Q x (0, oo)) a very weak ancient solution to equations 

(HIDO if 

/ u(x,t)(d t 4> + A(f))(x,t)dxdt = (1.4) 

-oo J Q 

for any ^ 6 C C M (0 x (-oo, 0)) with div = 0, 0|enx(-oo,o) = 0; and 

[ u(x,t)Vip(x,t)dxdt = (1.5) 
-oo Jn 

for any G C™{U x (-oo, 0)). 

Solutions defined in this way are often called very weak solutions in the literature and we 
also use this terminology. For smooth solutions the definition coincides with the usual one, 
as one can easily check by integration by parts. Our main result is as follows: 
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Theorem 1.1. Let u be a bounded very weak ancient solution to equations U.l\) . fT!fy . in 
the sense that u satisfies equations ( [1.^ ,( 17731) . Then u is given by ( II. 3\) . 

Remarks: The results are essentially sharp. This is obvious in the cases when Q is R n or a 
bounded domain. In the case of a half space, one can take u = (ui(t, x n ), . . . , u n _i(t, x n ), 0), 
Ui verifies: d t Ui — d\ui = fi(t), Ui\ Xn=0 = for 1 < i < n — 1, where fi £ C£°(— oo, 0). Then 
u is a solution to equations fll.ip . fll.2p . This is the example given in [8]. In exterior domain, 
the decay rate we obtain is as good as that of the fundamental solution of steady Stokes 
equations. 

Our work is motivated by boundary regularity for the Navier-Stokes equations. The main 
interest is in the case of the half-space, the other case are included for completeness. The 
connection between regularity and Liouville-type theorems is of course classical. In the 
context of the Navier-Stokes equations it is discussed for example in [01 Ej. In a recent 
note [S] a bounded shear flow for unsteady Stokes equations is constructed which is not fully 
regular although the boundary value is zero. This example simplifies earlier constructions 
of [5] • The lack of boundary regularity in the time-dependent case is in contrast with the 
case of steady Stokes equations, see e.g. [I]. In the time-dependent Stokes equations and 
Navier-Stokes equations, one usually treats pressure as an auxiliary variable, determined by 
u. Such treatment is valid as long as we have some decay of u at spatial infinity. On the 
other hand, it has been known that in unbounded domains, if we do not assume decay of u, 
the pressure may act as an external force 'driving' the fluid motion, as in the case of [5]. In 
such situations, we lose boundary regularity even with the vanishing boundary conditions. 
In this context, our result could be understood as showing that the solutions in [5] are in 
some sense the only obstacle to full boundary regularity (in suitable solution classes). 

The paper is organized as follows. In section 2 we introduce some technical lemmas to be 
used below. Section 3 deals with the simple cases when Q is a bounded domain or the whole 
space. Section 4 and 5 deal with the more subtle cases when Q is a half space or an exterior 
domain. For the exterior domains, we use a standard extension argument together with some 
estimates of linear Stokes system. For the half space, which is the most interesting case, we 
use Fourier transform. There is also a proof based on duality arguments, which requires 
some additional point-wise estimates of solutions to linear Stokes system in half space. The 
estimates may be of independent interest, but the calculations are somewhat lengthy. This 
alternative proof will appear elsewhere. 

Notation. We will use standard notations. For example, Q will be one of the four types of 
domain in R n mentioned above. B r (xo) := {x G R n \\x — xq\ < r}, Q r (xo,to) := B r (x ) x 
(to — r 2 , t ), Q r '■= Qr(0, 0). C denotes an absolute positive number. C(a, X, ■ • ■ ) denote a 
positive number depending only on a, A, • • ■ , A denotes the closure of A, A d O means the 
closure of A is a compact subset of O, d{ := 

2. Some technical lemmas 

In the sequel we will make use standard mollifications. For completeness we include the 
following standard lemma: 
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Lemma 2.1. Let Q be as above, u e L™ t (Q x (— oo,0)). Take a standard smooth cutoff 
function rj(t) with supp t] <<= (0, 1) and J r] = 1. For each e > 0, we define u e as a distribution 
in Q x (—oo,0) in the following way, 

r r° i s _ t 

I u(x,t) / -rj( )(p(x, s)dsdxdt (2.1) 

-oo Jn J-od e e 

for any smooth <p with supp <fi d Q x (— oo, 0). 

Then u e is a bounded function with bounded distributional derivatives d^u € , k = 0, 1, 2, • • • . 
Moreover, we have the following estimates: 

\\d t it 6 1| (fix (-00,0)) < Ce~ fc ||^||L°°(fix(-oo,o))- (2-2) 

Proof and Remarks: The proof follows immeditely from well-known properties of convo- 
lution. We note that due to our special choice of the support of rj, the mollified function u e 
is still defined in Q x (— oo,0). It is clear from definition that u e converge weakly* to u in 
L°°(Q x (—00,0)). It is also clear that, after possibly changing the value of u e on a set of 
measure zero, the map t — > u e (-,t) is continuous from (— oo,0) to L P (K) for any K d= Q, 
1 < p < oo. 

Let u be a bounded distributional solution to the linear Stokes equations (II. ip in Qi with 
some distribution p. It is well known that we have regularity of u in x, for almost every t 
in Qi/2- We can not, however, expect to have any regularity in t for u, or any reasonable 
estimate on p in general, assuming only that u is bounded in Q\. This point is usually 
illustrated with the example where u(t,x) = f(t), p = — fit) ■ x. Here f(t) is bounded, but 
f'(t) can be arbitrarily large. On the other hand, if we assume some estimate on d t u, then 
we can improve estimates on p. The following lemma summarizes the above discussion. 

Lemma 2.2. Let u be a bounded distributional solution to linear unsteady Stokes equations in 
Qi. Let IMUj^Q!) < 1- Then for any multi-index a with \a\ > 0, \\dxU\\L^ t (Q 1/2 ) < C(n,a). 
If in addition \\d t u\\Loo (q 1 ) < M, then there exists a pressure field p(x,t) such that U.l\) is 
satisfied and \\d%p\\L™ t (Q 1/2 ) — C{a,M,n) for any multi-index a with \a\ > 0. 

Proof: For the first part of the lemma, note that the vorticity Uij := diUj — djUi, 1 < i, j < n, 
satisfies heat equation dtuiij — Auij — in Q\. Thus is smooth with all derivatives 
bounded by constants depending only on n in Q3/4. From the divergence free condition, we 
get Aui = — Y21=i Then the first part of the lemma follows from interior estimate of 

Laplace equations. For the second part, note that || Vp||£<x> t (Q 3/4 ) < C(n, M) from the as- 
sumption on d t u and first part of the lemma. Since we also have Ap = 0, the estimate follows. 



Remark: The pressure is only determined up to an arbitrary function of t. (If we change 
p to p + c(t), equation ( II. ip is not affected.) In estimates below we will usually assume a 
suitable choice of c{t). 

We shall need the following extension result (which is interesting in its own right) below. 
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Lemma 2.3. (Extension of divergence- free vector field) For any smooth compactly supported 
vector field g = (g±, ■ • ■ ,g n -\,0) in R n , there exists a smooth divergence free vector field 
4> = (0i, ■ ■ • , 4> n ) with supp <p m R\ such that (j)\ Xn =Q = and ^-\ Xn=0 = g. 

Proof: We seek (j) in the form of <f>i = YTj=i ^i w ih with some Wy G C%°(R+) and Wij = —Wji 
for 1 < i,j < n. Note that under such conditions on Wij, div = is automatically satisfied. 
To satisfy boundary conditions for </>, we need: 

^^L=o = forl<z<n, (2.3) 



dx a 
3 =i j 



and 



^a^ L=0 = ^ forl <*<^n=o. (2-4) 

j=l n J 

It is easy to verify that the n-th equation in (12. 4 p is also automatically satisfied once the 
rest of the equations in the above are satisfied. To satisfy equations (12.31) (12. 4p . we first 
require Wij\ Xn=0 = 0, for all 1 < i,j < n. Then equations (12.31) become ^^\ Xn=0 = for 
1 < i < n. We further require that Wij — if 1 < i, j < n — 1, then equations (12 .4p reduce 
to d Q2^ n \ Xn =o — 9i for 1 < * < 71 — 1. Summarizing the above analysis, it is sufficient to find 
w in £ C^°(i?") for 1 < i < n — 1, w in = —w ni satisfying 

= 

: \ for 1 < i < n - 1. (2.5) 





X n = 


dw in 


x„=0 


dx n 


d Win 


x„=0 


d 2 x n 



■ 9% 

It is clear that we can always find such Wi n . Thus (j) satisfying conditions in the lemma exists. 

We collect some facts about the operator |V| which will be used in our proofs. 

For / G S(R n ), we define |V|/(x) = (|^|/(0) V ( X ) where we have used Fourier transform 
f(0 = j-^n J Rn e~ lx< f(x)dx and inverse Fourier transform f(x) = -^72 J Rn e tx< f(^)dC 

One can write |V|/ = Y^j=i ~w~\®if = Sj=i Rjfyf: where Rj denotes the Riesz transform. 
Clearly |V| can be considered as a continuous operator from S(R n ) — > L l (R n ). By duality, 
we can extend |V| to act on L°°(R n ) according to the usual fromula (|V|/, </>) = (/, |V|0) 
for any G S(R n ). 

We recall the following obvious continuity result. 

Lemma 2.4. Let |V|: L°°(R n ) 1 — > S'(R n ) be defined as above. If u m G L°° , with u m 
converges weakly* to u in L°° (viewed as the dual of I/ 1 (i? n ) y ), then |V|it m converges to |V|u 
in S'(R n ) 

Proof: This follows directly from the definitions. 
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Recall the definition of Holder norm in R n : H^llo^CR") := YlT=o su P|^|<m su Pxe-R n \d^u{x)\-\- 
sup^i^ sup z ygfini ^.^y for any m > 0, < a < 1. The Holder space C m,a (R n ) 

is consisted of all u with ||u||c , ™> a < oo. We will use the following estimate: 

Lemma 2.5. |V| : C m+1 > a (R n ) \-> C m ' a (R n ) is bounded, form>l, < a < 1. 



Proof: This follows from the representation |V| = Ylj ^j^j an d the Schauder estimates for 
the Riesz transform. 



We will denote by |V'| the analogue of |V| acting only on the variables x\, . . . ,x n -\. For 
Schwartz function /, | V'|/(x', x n ) = (2n\^'\f(^',x n )) w (x'), where the Fourier transform and 
inverse Fourier transform are both with respect to the first n — 1 variables. From defi- 
nition, it is clear if f(x',x n ,t) e L 00 ^™" 1 x (x 1 ,x 2 ) x (tx,t 2 )), then |V'|/ G V'{R n - x x 
(x u x 2 ) x (t l5 1 2 )). Moreover, if d l x ,d*d™f G L°°, then |V'|<%,d£d™/ = d l x ,d*d™\V'\f in 
V'iR"- 1 x ( Xl ,x 2 ) x (t x ,t 2 )). 

For bounded harmonic functions in the upper half spaces, we have the following result 
(see also [12], for example). 

Lemma 2.6. Let f be a bounded harmonic function in the upper half space i?" ; we have 
(dnf + I V'|/)(x) = in the sense of distributions in i?" . 



Remarks: By classical regularity for harmonic functions and lemma 12.51 we see both d n f 
and | V'|/ are smooth functions in the interior of i?". 

Proof: Let P(x,y) be the Poisson kernel. By classical representation results there exists 
a g e L°°(/? n_1 ), such that f(x) = f R „~i P(x,y)g(y) dy. By approximation and continuity 
properties of |V'| we can assume without loss of generality that g is smooth and compactly 
supported. Applying Fourier transform in the x\, . . . ,x n -i variables, we have /(^',x„) = 
g(£')e - ^ ' Xn and the result follows. 

3. The cases VI = R n or a bounded domain 

In this section, we first deal with the (easy) cases when Q = R n or Q is a bounded domain. 
Recall that our goal is to show that bounded very weak ancient solutions to fll.ipfll.2l) are 
given by (1 1 . 3 [) . 



i. n = R n . 

In this case, it is not difficult to see that equations (ll.4IUIl.5p are equivalent to 



d t u — Au + Vq = 
div u = 




in R n x (-oo,0) 



in the sense of distributions for some q € V(R n x (— oo, 0)). 

For 1 < i,j < n, let LUy = djUi — diUj, then clearly d t uiij — Aujij = in R n x (— co,0). 
Since Uy are bounded in some negative Sobolev space, we immediatly get ojy are bounded 
functions from parabolic regularity. Thus are so called bounded ancient solution to 
heat equation, and consequently oOy =constants c^. Since u is divergence free, we get 



(i 



Attj = — Y^j=i dj u ij = in i? n x (—00, 0). Therefore u(t, x) = f(t) for some bounded mea- 
surable / a.e t. This completes the proof when Q = R n . 



2. f2 is a bounded domain. 

In this case our goal is to show that bounded very weak ancient solutions u to (II. ip (11. 2p are 
identically 0. We use a duality argument as follows. For any / £ C^°(fl x (0,+oo)), let <fi 
solve 



d^-Afi + Vq =f 
div 6 =0 



in Q x (0, 00), 



4>(;t)\an = 0. 

The existence, uniqueness and regularity of such solutions are well known, one can see 
e.g [2]. Moreover, we have lim^oo \\(f>(-, i)||L 2 (n) = (the decay is actually exponential). 
Take a standard smooth cutoff function 77 (t) with rj(t) = for t > 2. For any R > 0, let 
4>n(x,t) = r](—j^)(f)(x, —t) for t £ (— 00, 0). Then from equations (jl.4p(ll.5p we obtain 



= / / u(x, t)(d t (j)R + A(p R )dxdt 
J -00 Jn 







00 J Q 




u(x, t)(—d t (f) + A0)(ar, —t)r]( — ^-)dxdt 

R 



■ u(x,t)rj'(—-)(j)(x,-t)dxdt 

it ./_oo 7f7 - K 

u(x,t)f(x,—t)rj(——)dxdt — — / / u(x,t)rj'(——)(j)(x,—t)dxdt 
00 Jn R RJ-ocJn R 



+ / / u(x,t)r)(—)Vq(x)dxdt 
J -00 Jn R 



Using the fact that / is compactly supported in t, q is smooth in x, u is bounded and 
lim^oo ||<^(-, t) ||ii(n) = (since Q is bounded), we can send R — > 00 and obtain 

-0 f 

u(x, t)f(x, —t)dxdt = 0. 



-00 Jn 

Since / is arbitrary, we must have u = 0. 

4. The case Vt = R^ 

Now let us deal with the more subtle case when Q is a half space. In fact one can still 
use the idea of duality as in the case of bounded domains. In this case, however, one has 
to study the decay property of solution to the linear Stokes equations quite carefully. One 
also has to appropriately localize <p (assuming notations from the last section) since in (11.41) 
the test function <fi is required to be of compact support. The authors have obtained a proof 
using such a method, which will appear elsewhere. 

Here we take a different approach based on the Fourier transform in which the calculations 



7 



are simpler. 

Let u be as above, take a smooth mollifier rj(x', t) with supp rj d -Bi(O) x (0, 1) C x R t 
and J i] = 1. We define the mollified vector field u e similar as before, again by duality: for 
any smooth (j) with supp <fi <s R+ X (— oo,0), 

(u e ,(j)): = / / u(:r,i) / / e~ n r](- — , - — -)(f)(y',x n ,s)dy'dsdxdt. 

J-oo J R™ Jfi"" 1 J-oo 6 e 

Again similar as before, one can show u e is bounded with bounded distributional derivatives 
\d^d^,u e \ < C(k, a, n)e~ k ~\ a \\\u\\L<x> . We have the following result: 

Lemma 4.1. Let u be a bounded very weak ancient solution to ( ti.il ) ( TT^) in i?" x (— oo,0) ; 
u e is defined as above. Then u e is smooth with all derivatives bounded in R\ x (— oo,0). 
Moreover, u e still satisfies equations ( |i.^| ) ; ( TO]) and u e (-,t)\ Xn=0 = 0. 

Proof: From equations fll.4p f ll.5j) and definition of u e we see 

-o 



u e (x,t)(d t (j) + A<j>)(x,t)dxdt = (4.1) 

oo J n 

for any <fi G C^°(Q x (— oo, 0)) with div = 0, 0|enx(-oo,o) = 0; and 

/ ( u t (x : t)V^(x,t)dxdt = V (4.2) 
J — oo J n 

for any if) G C^{VL x (— oo,0)). These clearly imply 

d ' U '~ AU ' + Zu< lo} inf'Wx(-oo,0)). (4.3) 

We first show u e is smooth up to boundary {x n = 0}. From the differentiablity property of 
u e in x',t, we see q is well defined for each t G (oo, 0) modulo some c(t). Moreover, from 
Aq = and elliptic estimates, we know q is smooth in x away from boundary {x n = 0}. 
Now let us rewrite the n-th equation of (14.31) as 

d 

= d t u £ n - Au e n + d n q = -—(q - d n u e n ) - A x ,u e n + d t u e n . 

ox n 

Note that d n u e n = — J2*i=i ®i u i ls bounded up to x n = 0. We see q — d n u e n is bounded up 
to boundary, thus q is bounded up to boundary {x n = 0}. The same argument also shows 
V",q is bounded up to boundary. Use Aq = we obtain that q is smooth in spatial variables 
up to x n = 0. Then from d^u e = d t u e — A x iu e + Vq we get u G C 2 . By differentiating the 
equations in x n and applying similar arguments we obtain smoothness of u e . Next we show 
u e \ Xn= Q = 0. Since u e is smooth in i?", we can use equations (11.41) (11.51) and integration by 
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parts in equations (I4.ipfl4.2p to obtain: 



fj 

J-oo J R n 

f J 

J-oo J R 



u^ipdxdt 



u e (j)dxdt 



0. 



0. 



(4.4) 



(4.5) 



Clearly ip can be arbitrary smooth compactly supported function, thus w^| Xn= o = 0. By 
lemma |2~3| (fi\ Xn =o can be any smooth compactly supported vector field with zero n-th com- 
ponent, thus u\\ Xn= o = for 1 < i < n — 1. Therefore, u e \ Xn = 0. 

Now we can prove our main theorem in this section. 

Theorem 4.1. Let u be a bounded very weak ancient solution to equations / li.ij) / TO|) in 
R+ x (—oo,0), then we must have u(x,t) = u(x n ,t) and u n = 0. 

Proof: By the above results, it is clear we only need to prove our theorem in the case u(t, x) 
is smooth up to boundary, with all derivatives bounded and u\ Xn= o = 0. Then we see J^- 

is bounded, for 1 < i < n. Since Aj^- = 0, from lemma [2761 we get (d n + |V'|)^- = 0. 
Applying operator (<9 n + |V'|) to the n-th equation of f ll.il) . noting the commutativity of 
various Fourier multipliers (since u is smooth), we infer that (d n + IV'Dw,,, satisfies the 
heat equation in R n x (— oo, 0) with (d n + \V'\)u n bounded and (d n + | V / |)ti n | x „ = o = (since 
d n u n = — Y^=i ®i u i anc ^ u \x„=o — 0). Thus by Liouville's theorem for heat equation in a half 
space, we get (d n + \V'\)u n = 0, and consequently Au n = 0. Since we also have u n \ Xn= o = 0, 
we see u n = 0. Therefore = 0. Thus |V'|Jj- = for 1 < i < n. Again applying operator 
| V'| to the first n — 1 equations of f ll.ll) . we get \V\u' satisfies heat equation in R n x (— oo, 0) 
with (|V'|M')|a; n =o = and |V'|u' bounded. Using Liouville's theorem for heat equation in a 
half space again, we obtain | V'\u' = 0. Thus u'(t, x) = u'(t, x n ). Summarizing the above, we 
obtain u(t,x) = u(t,x n ), and u n (t,x) = 0. 



Let u be a bounded very weak ancient solution to f ll.ipfll.2p in an exterior domain Q (i.e, 
the complement of Q is homeomorphic to a ball), we show u(x,t) = f(t) + 0( oi=s ) with 
some bounded / and n > 3, in this section. More precisely we have the following theorem: 

Theorem 5.1. Let u be a bounded very weak ancient solution to equations ( fi.Jj) ( fI7^j) in 

VL x (—00,0), where Q C R n (n > 3) is an exterior domain with Q c C B R for some R > 0. 
|l u IUj° t — 1- Then there exists a function a(t) with \a(t)\ < 1 for almost every t such that 



For such purpose, we first mollify u in t variable as in lemma 12. 1[ it is clear that u € thus 
obtained still satisfies equations fll.4p f ll.5p . Our first goal is to show that u e is smooth in 
x (-oo,0) and u € \ d n,t<o = 0. 

Lemma 5.1. Let u and u 6 be as above. Then u e is smooth infix (— oo, 0) and w e |an ) t<o = 0. 



5. The case is an exterior domain 



u(x, t) — a{t)\ < 



C(n,R) 



for almost every \x\> AR and t < 0. 



(5.1) 
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Proof: Clearly u e verifies 

d t u e - Au e + Vg = 1 . . . . . 

dw =Q | infix (-oo,0). (5.2) 

In equations fll.4ftfll.5p . take test functions as r/(t)<j)(x), r](t)tp(x) respectively for smooth 
0, ijj with supp 0, supp i(j fl, div = 0, 0|an = and 77 6 C£°( — 00, 0). We obtain by 
integration by parts (and definition of u e ): 

I I {-d t u e <p + u € A<j))ri{t)dxdt = 0, 

J — 00 J f2 

/ I u e ■ Vtpr](t)dxdt = 0. 

J —00 J Q 

Since 77 is arbitrary, we get for any t G (—00, 0), 

/ -d t u e <f) + u € A<pdxdt = 0, (5.3) 
Jn 

u e -V^dxdt = 0. (5.4) 



n 

Take R > sufficiently large such that f2 c C Br(0). For fixed t < 0, we can find v G 
^(flnBa) satisfying 

-At; + Vp =-d t u e (-,t)\ 

divv =0 / m " RjBi? ' (5 ' 5) 

w|sn = 0, v\ 9Br = u e (-,t)\ 9BR . (5.6) 

Note in the interior of Q, u e is smooth by lemma l2~2l and definition of u e . The existence of v 
follows from well-known results of steady Stokes system, we only remark here that the usual 
no outflow condition required by existence theory is satisfied in our situation and can be eas- 
ily seen by setting -0 to be 1 in a neighborhood of dQ in equation (I5.4p . Set w = u e (-, t) — v, 
we claim w = in Q D B R . To prove the claim, take any <fi G C°°(f2 H Br) with div = 
and <t>\a(nnB R ) = 0, ip G C°°(f2 H -B R ), we write = 0i + 2 , "0 = "01 + ^2 with the following 
properties: div 0i = div 2 = 0, 0i, 2 , -0i, -0 2 are smooth; 1; -0x equal and ip in a 
neighborhood of <9f2, vanishes in a neighborhood of OBr respectively. The existence of such 
decomposition of i\) is clear. To obtain such this decomposition for 0, one can localize by a 
standard cutoff function vanishing in a neighborhood of OBr, then use Bogovski's theorem 
to deal with the divergence free condition, we omit the details here. With these decompo- 
sitions, equations fl5.3pfl5.4p . the definitions of v and the fact that u e is smooth away from 
dfl, we easily obtain: f n wAcj)dx = and J n wVifidx = 0. Thus by result in section 3, this 
implies w — 0. Therefore, u e (-,t) G C 1,1 I 2 (Br fl fi) and u £ \qq = 0. A simple boostraping 
argument gives smoothness of u € . The lemma is proved. 



Proof of main result of this section 

Let us first summarize the above results as follows: 
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u e is in C°°(Q x (— oo,0)) with all derivatives bounded (with bounds depending on e) and, 
u e satisfies 

d t u e - Au e + Vq = \ . . m , . 

divtf =0 j mnx(-oo,0), (5-7) 

u e (;t)\ dn = 0. (5.8) 

We extend u e to f? n by setting u e = in f2 c . It is not hard to see the extended u e satisfies 

d t u e - Au e + Vq = n 1 . Dn , „s /c „ A 

divn =0 / m ^ x (-°°.°)- ( 5 - 9 ) 

for the measure [i = f e (x, t)da, where f e = — qn on d£l and da is the surface measure of 
dQ. We set 

x-y 



v*(x,t)= I Pe A{t ~ s) n(-,s)ds = I -j-J k(-=L)ny,s)da(y)ds 

A 1 ~ s ) ' Jan v t_s 



where P is the Helmholtz projection to divergence free vector field and fc(-) is the kernel of 

C(n) 

(i+\y\) n ' 



Pe A . Thus \k(y)\ < frrfb ■ Simple calculations show v e verifies the following estimates: 



\\v e (;t)\\LHB 2R ) < C(R,n,e), and \v c (x,t)\ < |W(M)| < for |x| > 2R. 

In the above calculations n > 3 is important. One can check if n — 2 the integral might 
diverge. Clearly w e := u e — v e satisfies 

d t w e -Aw e + Vq — 1 . Dn , nN /c m\ 

div J =0 ) in ^x(-oc,0). (5-10) 

Thus u> £ = a e (t) and consequently u e = v e + a e (t). At this stage, we would like to pass e 
to zero. The decay estimate for v € , however, depends on e (since the bounds of f e depends 
on e). Thus we must first remove this dependence. To do this, let us consider vorticity 
ulj = diUj—djUj for 1 < i, j < n. ufj satisfy dtufj—Aufj = in (R n \B 2 R) x (— oo, 0). By inte- 



nt: scaling nivaiifciiiut; uj^yx, l) — t uu^yivix, ivl 2 

and the L°° bound on u, we easily conclude: 



rior regularity of solution to heat equation, the scaling invariance wL(x, t) — > w-^Mx, M x) 



u>lj is smooth in (R n \B 3 ji) x (— oo,0) with \d t <jJ^{x^t)\ < C Zl? , \x\ > 3R, the estimate being 
independent of e. 

Now fix e = eo > 0. From estimates of u e ° above, we know \u^°(x,t)\ < ^fe^ff , \x\ > 3R. 
By estimates of d t u>ij and definition of ui^, we see \uij(x, t) — uj^(x, t)\ < C ^° 3 ' ir> , |x| > 3R. 
Therefore, \uij(x,t)\ < C(n, eo, R)(t^ + up), \x\ > 3R. In fact, by a boostrap argument 
(better decay estimate of u)y improves estimate of d t Uij which in turn implies better decay 
estimate of tOij), one can upgrade this estimate to \u>ij(x, t)\ < ^prp ■ Thus \u e (x,t)\ < 
(n,eo,R) ^ ^ now independent of e) for \x\ > 3R. For each fixed t G (— 00, 0), u\ solves 



Aut = -J2 d Mj ™R n \B 3R . (5.11) 
j'=i 
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From this and the boundedness of u, it is not hard to see u\ verifies the following bound: 
\u e (x,t) — a e (t)\ < C '|"|'^°'f ? ' > for |x| > 4f?. Passing e — > the conclusion of this section is 
reached. 
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